Measurement of the Strong Coupling Constant 
from Inclusive Jet Production at the Tevatron pp Collider 



We report a measurement of the strong coupling constant, as{Mz), extracted from inclusive jet 
production in pp collisions at ^/s =1800 GeV. The QCD prediction for the evolution of as with jet 
transverse energy Et is tested over the range 40 < Et < 450 GeV using Et for the renormalization 
scale. The data show good agreement with QCD in the region below 250 GeV. In the text we 
discuss the data-theory comparison in the region from 250 to 450 GeV. The value of as at the mass 
of the Z" boson averaged over the range 40 < < 250 GeV is found to be as{Mz) = 0.1178 ± 
0.0001 (stat)^y pggg(exp. sy st ) . The associated theoretical uncertainties are mainly due to the choice 

of renormalization scale {t.^J a-nd input parton distribution functions (5%). 



PACS numbers: 12.38.Qk, 13.87.Ce 
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Jet production at hadron colliders provides an excel- 
lent opportunity for testing the theory of strong interac- 
tions, Quantum Chromo dynamics (QCD) [Q. QCD has 
achieved remarkable success in describing hadron inter- 
actions at short distances (large momentum transfers), 
owing to the property of asymptotic freedom . Asymp- 
totic freedom predicts a logarithmic decrease of the cou- 
pling strength, as{^J')^ as the momentum scale /i char- 
acterizing a process increases. Processes with large mo- 
mentum transfer can then be described by an expansion 
in powers of Q;s(/Lt). The value of as, a free parameter of 
QCD, is one of the fundamental constants of nature. Its 
determination is the essential measurement of QCD, and 
the observation of its evolution, or running, with momen- 
tum transfer is one of the key tests of the theory. At e+e^ 
colliders as has been measured from the fragmentation 
functions Q, event shapes [Q, jet production rates ||] 
and T lepton decay properties Q. In lepton-hadron col- 
lisions, as has been measured from scaling violations [Q, 
jet production rates ||] and momentum sum rules |^ . A 
precise value for as has also been obtained from a global 
fit to properties of the Z° boson measured at the CERN 
LEP and the SLAC SLC e+e" colhders and the W bo- 
son and top quark masses jlOj. A review of these and 
other as measurements can be found in . In this let- 
ter, we present a measurement of as from the inclusive 
jet cross section in pp collisions over the jet transverse 
energy {Et) range from 40 to 450 GeV. 

This measurement is based on a data sample of inte- 
grated luminosity 87 pb~^ collected by the Collider De- 
tector at Fermilab (CDF) during the 1994-95 run (Run 
lb) of the Fermilab Tevatron pp collider operating at 
= 1.8 TeV. The CDF detector is described else- 
where |l^. Details of the measurement of the inclusive 
jet differential cross section can be found in jl^. Briefly, 
jets are reconstructed using an iterative fixed cone algo- 
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rithm with a radius R = (Ary^ + Acf)'^) — 0.7, where 
7] = — ln(tan|) is the pseudorapidity, evaluated from the 
angle 6 between the centerline of the jet cone and the 
proton beam line, and is the azimuthal angle. The 
inclusive jet cross section includes all jets in an event 
in the pseudorapidity range 0.1 < \ri\ < 0.7. The mea- 
sured spectrum is corrected for the calorimeter response, 
resolution and the underlying event energy using an iter- 



ative unsmearing procedure which changes both the en- 
ergy scale and the normalization simultaneously. The 
value of as is determined by comparing the jet cross sec- 
tion with the next to leading order (NLO) perturbative 
QCD predictions |l^. In the Et region studied, the non- 
perturbative contributions to the inclusive jet cross sec- 
tion are estimated to be negligible ||l^. The procedure 
of extracting as can be summarized by the equation 



da 
dE^ 



= al{nR)X''°\^F,ET) [1 + as{iJLR)ki(p.Fi, ^jip , Et)] 

(1) 



where -j^^ is the transverse energy distribution of 
the inclusive jets, iia and jip, related to Et by a 
scale factor, are the renormalization and factoriza- 
tion scales, al{iiB.)X'^^\fXF, Et) is the leading order 
(LO) prediction for the inclusive jet cross section, and 

al{p,Ti)X^'^\^F,ET)ki{^R,^F,ET) the NLO contribu- 
tion. Both X^^\iif,Et) and ki{^ii, plf,Et) are calcu- 
lated with the JETRAD Monte Carlo program |l^ based 
on the techniques described in ||l^ and the matrix ele- 
ments of |l^. NLO QCD predictions for the inclusive jet 
cross section are also available in ||l^, 20 1 and agree well 
with those of jetrad. All calculations are performed 
in the modified minimal subtraction, MS , scheme . 
The JETRAD Monte Carlo program generates events with 
weighting factors, so that the jet clustering algorithm and 
Et and rj cuts, mimicking the experimental requirements, 
are directly applied to the final state partons. The jet 
clustering in jetrad is governed by a cone radius R and 
a phenomenological parameter 7?.sep (default value =1.3), 
introduced to match the experimental efficiency of iden- 
tifying overlapping jetsJTst . If two partons are more than 
TZsep X R apart, they are identified as two distinct jets, 
otherwise they are merged into a single jet. 

The inclusive jet data are divided into 33 Et bins, 
from which we obtain statistically independent measure- 
ments of as for 33 different values of /i^. The as values 
derived ior = fj,F — Et using CTEQ4m ||2^ parton dis- 
tribution functions (PDFs) are presented in Fig. |l|. For 
Et <250 GeV, there is good agreement with QCD pre- 
dictions for the running of the coupling constant. The 
behavior of as at higher Et values is a direct reflection of 
the excess observed in the inclusive jet cross section [n3|. 
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FIG. 1: The strong coupling constant as a function of Et 
for fiR — Et measured using CTEQ4m parton distributions. 
The shaded area shows the experimental systematic uncer- 
tainties. The curved line represents the NLO QCD prediction 
for the evolution of asiEr) using as(Afz) = 0.1178, the av- 
erage value obtained in the region 40 < Et < 250 GeV. The 
as{Mz) extracted from Qs(i?T) is shown in the inset along 
with the weighted average as the horizontal line. 



The discrepancy with the NLO QCD predictions in this 
region, though not well understood, may be accommo- 
dated by the fiexibihty allowed by the world data in de- 
termining the high-x gluon component in the parton dis- 
tributions p2| . 

The measured values of as {(J-r) are evolved to the mass 
of the Z° boson, Mz, by using the solution to the 2-loop 
renormalization group equation 



as{Mz 



1 - as(/i_R.)(&o + hasi^iR)) ln{fj,R/Mz) 



(2) 



bn = 



33 - 2n 



f 



67r 



306 - 38n^ 
- 24^2 . (3) 



where f^/ is the number of active flavors, which is equal 
to five (six) for smaller (larger) than the top quark 
mass. The values of as (M^) for all 33 measurements are 
shown in the inset of Fig. m. Averaging over the range 
40-250 GeV, we obtain 

as{Mz) = 0.1178 ± O.OOOl(stat). 

Inclusion of the data with Et >250 GeV results in an 
increase of the average value by O.OOOl. 

The running of as is tested by verifying if as{Mz) 
is independent of the energy scale Et at which the jet 



10 Tligh p.,- Pions Response 
nr- 



-10r 



LU 10 Energy Scale 



OE 



5-1 Or 



10 HJnderlying Event 
" 



0) 
Q. 



io: 



10 H-iesolution 
-lOr 



1 Tow p.,. Pions Response 





-lOr 



10 INeutral Pion Response 
01- ~ 



-lOr 



10 ^Fragmentation 





-io: 



10 ^vJormalization 
nt "~ 

-lOr 



40 130 220 310 400 40 130 220 310 400 
Transverse Energy (GeV) 

FIG. 2: Experimental systematic uncertainties for Os mea- 
surement (the lines are 1 standard deviation contours), with 
CTEQ4m as input PDF and = jiF — Et- 



cross section is measured. The 27 values of as{Mz) ob- 
tained from the data in the jet Et range 40-250 GeV are 



fit to the linear function Po + A x {Et/E^ — 1) with 
E^^ = 92.8 GeV. The fit yields Pq = 0.1176 ± 0.0003 and 
Pi = 0.0003 ± 0.0003 with x^/d.o.f. = 1.3, showing that 
as{Mz) is independent of Et within one standard devia- 
tion. This non-trivial result demonstrates the correctness 
of the QCD prediction for the evolution of as over the 
above range. 

The experimental systematic uncertainties on the value 
of as{Mz) are derived from those on the inclusive jet 
cross section. For each source of systematic uncertainty 
described below, except normalization, the inclusive jet 
cross section was re-evaluated by varying the correspond- 
ing parameter in the detector response by 1 a. For the 
normalization uncertainty it was changed by a scale fac- 
tor jl^. These uncertainties were propagated to as{Mz) 
by repeating the procedure described above using the 
spectra given in Table VI of Ref. ||l^. The results are 
shown in Fig. The deviations of as{Mz) for each spec- 
trum from the central value are given in Table H 

The dominant experimental systematic uncertainty in 
the inclusive jet cross section measurement is due to the 
calorimeter response to jets. The detector response and 
jet energy corrections are derived from a combination 
of test-beam data and Monte-Carlo simulations. The 
calorimeter response to charged pions was evaluated sep- 
arately for high and low transverse momentum (Pt) pi- 
ons from test beam data and isolated charged tracks from 
pp data with an uncertainty of ±5% for Pt < 5 GeV, 
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TABLE L Experimental systematic uncertainties on Qs(Mz) 
extracted using CTEQ4m parton distribution functions. 



Source of Uncertainty 




^ as ^ 




16 


Calorimeter high Pt pion response 


+ U.UUOO 
-0.0055 


+ 3.1 
-4.7 




0.14 


Calorimeter low Pt pion response 


+0.0027 
-0.0033 


+2.3 
-2.8 






Energy scale 


+0.0030 
-0.0030 


+2.6 
-2.6 


CO 

a 


0.12 


Neutral pion response 


+0.0016 
-0.0021 


+ 1.4 
-1.8 


c 
re 




Underlying event energy 


+0.0025 
-0.0027 


+2.1 
-2.3 


U) 

c 


0.10 


Jet fragmentation functions 


+0.0046 
-0.0044 


+3.9 
-3.7 


o 
O 




Jet energy resolution 


+0.0015 
-0.0017 


+ 1.3 
-1.4 


D) 
C 




Normalization 


+0.0022 
-0.0023 


+2.0 
-1.9 


■q. 





±3% for 5 GeV<PT<15 GeV and t^-S^ ^j. > 
GeV. During the run, the calorimeter response was mon- 
itored by using muons, isolated particles and the mea- 
sured inclusive jet cross section. The response was found 
to be stable within ±1%. The electromagnetic calorime- 
ter was calibrated using electrons from pp interactions. 
The associated uncertainty, labeled in Fig. || as neutral 
pion response, arises from the modeling of calorimeter 
response to very low energy electrons. The underlying 
event energy (non-jet energy contribution to the jet Et) 
was measured from minimum bias data and its mean 
value was varied by ±30% to evaluate the effect on the jet 
cross section. The error from the jet fragmentation func- 
tions is due to the extrapolation of the track momentum 
and multiplicity distribution to the high Et region and 
from uncertainties in the track reconstruction efficiency. 
The detector jet energy response has a Gaussian shape 
with exponential tails on both the high and low sides and 
a resolution with an uncertainty of ±10%. Finally, the 
overall normalization of the inclusive jet cross section has 
an uncertainty of ±4.5%, dominated by the contribution 
from the total cross section measurement. Summing in 
quadrature all the above uncertainties after propagation 
to as{Mz) yields a total experimental systematic uncer- 
tainty of±0:88^i. 

The theoretical uncertainties are mainly due to the 
choice of renormalization and factorization scales and 
parton distribution functions. The scales fxp and fiR 
are expected to be of the same order as the character- 
istic scale of the process, which in this case is the jet 
Et- We have evaluated the changes in as{Mz) result- 
ing from independently varying and from Et/^ 
to 2Et and found that the largest changes occur for 
Mfl. — I^F- For all results presented in this letter the 
two scales were set equal. The sensitivity of the mea- 
sured as{Mz) to changes in these scales is indicated by 
the shaded band in Fig. |^(a). Over the Et range from 
40 to 250 GeV, the shift in as{Mz) induced by changing 
the scale from Et/2 {as{Mz) = 0.1129 ± 0.0001) H to 
2Et {as{Mz) = 0.1249 ±0.0001) is approximately +^^^, 
independently of Et- 

The coefficients Xg and fci in Eq.( |^) depend on the 
PDFs, which are obtained from global fits to deep inelas- 
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FIG. 3: Uncertainties in aa{Mz) due to the renormalization 
scale fi, (a), and parton distribution functions, (b). 



tic scattering (DIS), Drell-Yan production and other col- 
lider data, including the inclusive jet results from Tcva- 
tron. Each PDF set has an associated strong coupling 
constant, af^^- The gluon PDF {G{x)) determined 
from the fit is highly correlated with a^^^ because in 
equations describing all the processes used, the G{x) is 
always accompanied by af^^- To calculate the above 
coefficients, the PDFs are evolved using af^^. For this 
procedure of measuring to be valid, the extracted 
value of as should be consistent with the input af^^, 
although not necessarily equal. The variation in parton 
distributions, especially in the gluon distribution, allowed 
by the world data was studied by the CTEQ Collabora- 
tion by fixing the value of af^^ to 0.110, 0.113, 0.116, 
0.119 or 0.122, with resulting of 1388, 1323, 1323, 1388 
or 1543 for 1297 non-jet data points |2^. We use the 
CTEQ4a series to study the as{Mz) dependence on the 
PDFs. In addition, we have studied as{Mz) using PDF 
sets which do not include Tevatron jet results, such as the 
MRST(gt) set [|4j, the mrsa' series and two iviRS-R 
sets . The calculated by comparing the data with 
the theoretical prediction in the restricted range of 40- 
250 GeV, is used to quantify the agreement. The minimal 
X^/d.o-f- = 1.38 is obtained for CTEQ4m (af ^^=0.116), 
therefore we use this PDF in our final fit. Excluding the 
PDFs which have obvious disagreement (x^/d.o./. > 5), 
we estimate the uncertainty on the as{Mz) due to the 
PDF choice to be ±5%. 



Finally, the variation of TZsep, the jet clustering pa- 
rameter, from 1.3 to 2.0 results in a 5-7% normalization 
change of the inclusive jet cross section. The correspond- 
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ing variation in the as{Mz) measurement is 2-3%. 

To explore the flexibility in the gluon distribution 
at high Et: a special PDF set, CTEQ4hj, was gener- 
ated by including CDF jet data in the global fit with 
higher statistical weight assigned to high Et points and 
a new parameterization of the gluon distribution [^ . 
This PDF yields good agreement between Tevatron data 
and theoretical predictions. Using this set, we obtain 
as{Mz) = 0.1185 ± 0.0001, averaged over the entire Et 
range. 

In conclusion, we have tested the evolution of the 
strong coupling constant as using the inclusive jet cross 
section data from pp collisions at y/s — 1800 GeV. Our 
results demonstrate that for Et in the range of 40-250 
GeV with = Et the running of is in good agree- 
ment with QCD predictions. The value of as expressed 
at the Z*^ boson mass is found to be 

as{Mz) = 0.1178 ± 0.0001(stat)j:° °[J^^(exp.syst) 

This value is in good agreement with the world average 
asiMz) = 0.1181 ± 0.0020 The theoretical uncer- 



tainties associated with the choice of parton distribution 
functions (~ 5%) and the choice of the renormalization 
scale (l4%) are comparable to the experimental system- 
atic uncertainty. 
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